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We consider a simulated population of isolated Galactic neutron stars. The rotational frequency of 
each neutron star evolves through a combination of electromagnetic and gravitational wave emission. 
The magnetic field strength dictates the dipolar emission, and the ellipticity (a measure of a neutron 
star's deformation) dictates the gravitational wave emission. Through both analytic and numerical 
means, we assess the detectability of the Galactic neutron star population and bound the magnetic 
field strength and ellipticity parameter space of Galactic neutron stars with or without a direct 
gravitational wave detection. While our simulated population is primitive, this work establishes a 
framework by which future efforts can be conducted. 



I. BACKGROUND AND MOTIVATION 

Isolated neutron stars with non-axisymmetric deforma- 
tions will continuously emit gravitational waves as they 
rotate pQ. Neutron stars also have strong dipolar mag- 
netic fields that accelerate particles to relativistic ener- 
gies . Since these neutron stars can lose energy through 
the emission of electromagnetic and gravitational radi- 
ation, their rotational frequency slowly decreases over 
time. The gravitational wave strain amplitude of rotating 
neutron stars has a strong dependence on the star's rota- 
tional frequency. Though no gravitational wave detection 
has yet been reported, rapidly rotating isolated Galactic 
neutron stars are one of the most promising sources of 
continuous gravitational waves for ground based gravita- 
tional wave detectors. 

Attempts to assess the detectability of gravitational 
waves from the Galactic neutron star population began 
with rough analytic estimates. An argument presented 
in [3] by Thorne but credited to Blandford models the 
Galactic neutron star population as a uniformly pop- 
ulated two-dimensional disk of gravitars (neutron stars 
with gravitationally dominated frequency evolutions) all 
born at very high frequencies. Using this simplistic 
model, he estimated a rough upper bound on the possible 
gravitational wave strain amplitude from a Galactic neu- 
tron star, /i max ~ 10 -25 [3]. Blandford also surprisingly 
observed that the maximum gravitational wave ampli- 
tude is independent of the size of the star's deformation 
and rotational frequency. His argument was revised in 
[4] and again in [5], which both found /i max ~ 10 -24 . 

This work was followed by more comprehensive at- 
tempts to assess the detectability of the Galactic neu- 
tron star population through population synthesis. If 
the neutron star population can be accurately simulated, 
then the detectability of Galactic neutron stars can be 
determined. In [Bj Palomba was the first to assess the 
detectability of a simulated gravitar population by first 
and second generation gravitational wave detectors. He 
incorporated realistic spatial, age, birth frequency, and 



kick velocity distributions, as well as a possible elliptic- 
ity distribution (though this is still largely unconstrained 
[6]). He estimated the fraction of the neutron star popu- 
lation that would likely have to be gravitars in order for 
first or second generation detectors to make a direct grav- 
itational wave detection. Continued efforts by Knispel 
and Allen extended Blandford's argument to a simulated 
gravitar population similar to Palomba's [5] . They found 
that the maximum gravitational wave strain amplitude 
does have a strong dependence on the star's frequency 
and size of deformation when considering a more real- 
istic neutron star population. They set upper bounds, 
which depend on the population's ellipticity (a measure 
of a star's deformation) and rotational frequency, on the 
gravitational wave strain amplitude of the nearest source. 



In this paper, we include electromagnetic emission as 
well as gravitational wave emission in the frequency evo- 
lution of neutron stars and investigate its effect on the 
population's detectability. We use the simulated neutron 
star population in [5] and assign every neutron star a 
dipolar magnetic field as well as an ellipticity. We then 
allow each star's frequency to evolve through the emission 
of both gravitational and electromagnetic radiation. The 
paper is organized as follows. In Section [il] we review the 
spin and strain evolution of neutron stars and revisit the 
upper bounds from the gravitar case. In Section |IH| we 
outline a Monte Carlo simulation used to assess the de- 
tectability of the Galactic neutron star population. The 
results are then used to bound the magnetic field strength 
and ellipticity parameter space of isolated neutron stars 
with or without a direct gravitational wave detection. 



In Section IV we present a rough analytic argument to 
which we compare our numerical results. In Section [V] 
we summarize our main results. 
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II. SPIN AND STRAIN EVOLUTION OF 
NEUTRON STARS 



birth frequency vq 3> v to spin down to a frequency v) 
for gravitationally dominated emission is 



We use the simulated neutron star population from [5] 
to assess the detectability of gravitational waves emitted 
by isolated Galactic neutron stars. It is important to 
note that, while the simulated population does not ex- 
plicitly include recycled millisecond pulsars, it does not 
necessarily exclude them either. Each star in our popu- 
lation is assigned a birth frequency, initial position, kick 
velocity, and age. Stars are then independently evolved 
through the Galaxy's gravitational potential (see [5j for a 
more detailed description of the population). Therefore 
an old star that has been recently recycled can just be 
thought of as a young star born with a high frequency. 
We also consider a large enough range in magnetic field 
strength to accommodate recycled pulsars. In this sec- 
tion, we review methods to find the spin frequency and 
gravitational wave strain amplitude of each star in our 
simulated population in order to assess its detectability. 

If neutron stars only lose energy through gravitational 
and electromagnetic emission, their rotational frequency 
evolution is given by 
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in cgs units [7] 00. Here, G is the gravitational con- 
stant, c is the speed of light, v is the star's rotational 
frequency, R is the star's radius, I = kMR 2 is the mo- 
ment of inertia about its rotational axis with M being 
the star's mass and k w 2/5 [2 , e = (I± — is the 

ellipticity with I\ and 1% being the moments of inertia 
about the star's other two principle axes, B is the dipolar 
magnetic field strength at the magnetic equator, and a is 
the angle between its magnetic pole and its axis of rota- 
tion 1 . We choose the canonical values of R = 10 km and 
M = IAMq for all neutron stars [2J. Because we only 
concern ourselves with order of magnitude estimates, we 
set sin 2 a = 1 . 

Equation [T] can be solved analytically for v(t, in the 
limits where e = or B = 0. If B = 0, a neutron star 
will only emit gravitationally. Its frequency is 



(2) 



where t is the neutron star's age, v$ — v{t = 0) 
is the neutron star's birth frequency, and 7 gw = 
-5127r 4 Gic- 5 e75. Equation |] is a good approxima- 
tion for the frequency of a gravitar. The characteristic 
timescale (the approximate time for a neutron star with 



The second term in Equation[T] which is the frequency evolution 
due to electromagnetic emission, is derived from the simple model 
of a rotating dipole. In [9] Spitkovsky corrects this term such 
that a neutron star will still emit electromagnetically even if its 
magnetic pole and rotational axis are aligned. 
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If e = 0, a neutron star will only emit electromagnetically. 
Its frequency is 



u(t, vq) = (y 2 - 27di p £) 
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where 7di P = -8n 2 R 6 c~ 3 I~ 1 B 2 sin 2 a/3. Equation[i]is a 
good approximation for the frequency of a neutron star 
whose evolution is dominated by electromagnetic emis- 
sion and whose characteristic timescale is 
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While e is unknown, the dramatically different timescales 
between Equations [3] and [5] illustrate the difficultly in 
detecting isolated neutron stars: stars with reasonable 
magnetic fields spin down to low frequencies too rapidly 
to detect. Therefore, gravitational wave detectors will 
likely only detect neutron stars with small magnetic fields 
or young neutron stars that have not yet spun down to 
low frequencies. 

Not all neutron stars will have their frequency evolu- 
tion dominated by either gravitational or electromagnetic 
emission. For these stars, i> cannot be integrated over 
time to solve for an analytic solution for v(t, Vo\e, B). 
However, [S] shows that Equation [l] can instead be in- 
verted to solve for t{v, v^\e, B). Following [5], we rewrite 
Equation [T] as 



v = 7 gw i/ +7di P ^ 
= 7di P (7^ 5 + , 
where 7 = 7gw / 7 di p . Equation [7] can be solved for 
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If vq, t, 7gw , and 7 di p are known, Equation[H]can be solved 
numerically to find v using root finding techniques [10] . 

The strain amplitude of gravitational waves emitted by 
a neutron star at a radial distance r away from a detector 
is given by 
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assuming that the neutron star's sky location intersects 
a line normal to the plane of the detector arms and its 
axis of rotation is parallel to that line (optimal mutual 
orientation). Since we only concern ourselves with or- 
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der of magnitude estimates, we assume optimal mutual 
orientation for all neutron stars [5], which overestimates 
the detectable amplitude by about a factor of four on 
average. 

For a population of neutron stars whose radial distance 
from Earth r, age t, birth frequency vq, ellipticity e, and 
magnetic field strength B are known, Equation [2j |4j or 
[8] can be used to determine each star's spin frequency 
v. Equation [2] is used when 7 gw z /5 3> 7dip* y3 , which we 
conservatively choose to be when 7 > 40 s 2 ; Equation [4] 
is used when 7 gw t /5 7dip^ 3 , which we conservatively 
choose to be when 7 < 4 x 10~ 9 s 2 ; Equation [5] is used 
otherwise 2 . Equation [9] can further be used to determine 
each star's gravitational wave strain amplitude h as mea- 
sured in our detector. We compare each star's frequency 
and strain amplitude to a scaled gravitational wave de- 
tector's noise curve in order to assess the detectability of 
the neutron star population. We explain how we derive 
the scaling factor in Section [TTT] 

While Equation [9] for the gravitational wave strain am- 
plitude h does not depend explicitly on the magnetic 
field, B does help to determine v through Equation [T] 
There are two related effects. First, Figure [l] shows that, 
all other things being equal, neutron stars with large 
magnetic fields will spin down to low frequencies (high 
periods) much faster than neutron stars with small mag- 
netic fields. Consequently, large magnetic fields will re- 
sult in smaller and smaller gravitational wave amplitudes 
over time. Second, since gravitational wave detectors are 
sensitive to finite frequency ranges, neutron stars with 
large magnetic fields will rapidly spin through a detec- 
tor's sensitive frequencies, which makes them less likely 
to be detected. Therefore, neutron stars with small mag- 
netic fields are more likely to be detected than neutron 
stars with large magnetic fields. In this way, assuming 
we know the population's ellipticity, we can place lower 
bounds on the magnetic field of neutron stars in the ab- 
sence of a gravitational wave detection. 

We can gain intuition into the detectability of Galac- 
tic neutron stars by setting B = 0. This places an upper 
bound on h for fixed e values. In Figure[2j we plotted the 
maximum gravitational wave strain amplitude /i max ver- 
sus gravitational wave frequency / = 2u of the simulated 
neutron star population presented in [5] with B = and 
e = HT 9 , HT 8 , 1CT 7 , and 1CT 6 . A single point (/, h max ) 
corresponds to the population's maximum gravitational 
wave amplitude h max measured in the frequency band 
[/, ef] where e is Euler's number. 3 



2 To determine the two 7 cutoffs, we assume that one term will 
dominate over the other if it is at least three orders of magnitude 
greater than the other. Equation I2I can be used when 7» 1/v 2 . 
Therefore, we choose 7 > W^/u = 40 s 2 for v = 5 Hz. Equa- 
tion [4] can be used when 7 <C 1/u 2 . Therefore, we choose 
7 < lQ~ 3 /u 2 = 4 X 1CT 9 s 2 for v = 500 Hz. 

3 To find /imaxi we considered 200 logarithmically spaced overlap- 
ping frequency bands and constructed histograms for the strain 
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FIG. 1: This figure shows the period evolution of neutron stars. 
The dots represent observed pulsars from the ATNF catalog 
The thin, negatively sloped solid contours are lines of constant 
magnetic field strength (labels on the left), and the thin, positively 
sloped dotted contours are lines of constant characteristic age (la- 
bels on the right) assuming only electromagnetic emission. The 
thick, solid lines with square ticks track the period evolution of 
a neutron star that emits both electromagnetic and gravitational 
radiation. These lines, which are labelled by the logarithm of the 
star's magnetic field in units of Gauss, correspond to a neutron 
star with e = 10~ 7 (all lines) and B = 10 8 5 - fO 11 with steps 
of 1/2 dex. The square ticks represent logarithmic steps in age. 
The first tick labels t = 0, and the subsequent ticks range from 
t = 10 4 — 10 9 yrs. The thick, dashed lines, which are labelled 
by the logarithm of the star's magnetic field in units of Gauss, 
are characteristic aLIGO sensitivity curves for neutron stars with 
e = 10~ 7 located 100 pc away from Earth. Neutron stars below 
their associated aLIGO sensitivity curve are undetectable. Neu- 
tron stars with large magnetic fields spin down to low frequencies 
(high periods) much faster than stars with small magnetic fields; 
consequently, they spend less of their lives emitting gravitational 
waves with frequencies that aLIGO is most likely to detect. 

Our numerical result in Figure [2] is consistent with the 
result in [5], which was derived using a semi-analytical 
integration technique. Considering a distribution in fre- 
quency and a three-dimensional spatial distribution re- 
sults in a clear dependence of h max on both frequency 
and ellipticity [5]. The effect of the frequency distribu- 
tion manifests itself in the overall shape of the four curves 
in Figure [2] Since stars with large ellipticities spin down 
much faster than stars with small ellipticities (Equation 
[3]), a neutron star population with large ellipticities more 
densely populates low frequency bands than a neutron 
star population with small ellipticities. Therefore, while 
each curve in Figure [2] has a similar shape, the large el- 



amplitude from the neutron stars in each band. We then solved 
for hmax using a linear fit in logio-space to the tail (largest h 
values) of each histogram. We used this method to minimize 
statistical fluctuations. 
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FIG. 2: We plot the maximum gravitational wave strain ampli- 
tudes h ma x as a function of gravitational wave frequency / = 2v of 
a population of gravitars (neutron stars with B = 0) with e = 1Q~ 9 , 
1CT 8 , 10 -7 , and 10 -6 . A single point (/, h max) corresponds to the 
population's maximum gravitational wave amplitude h 

max mea- 
sured in the frequency band [/, ef] . We used the gamma initial 
radial distribution from [5] to simulate the neutron star popula- 
tion. 



lipticity curves arc shifted to the left relative to the small 
cllipticity curves. The subtle kink in the e = 10 -9 and 
10~ 8 curves between the nearly flat, high frequency re- 
gion and the more positively sloped, low frequency re- 
gion corresponds to a kink in the population's frequency 
distribution, which is described in [5]. 4 The e = 10~ 7 
and 10 -6 curves also exhibit the same behavior but at 
smaller frequencies than those plotted. If we had consid- 
ered a two-dimensional spatial distribution, /i max would 
have been independent of the ellipticity in the region to 
the right of the kink. Here, the frequency distribution is 
in a nearly steady state. Considering a three-dimensional 
spatial distribution breaks this degeneracy between ft max 
and the population's ellipticity. Note that the gravita- 
tional wave strain amplitude will decrease when magnetic 
fields are considered. 



III. NEUTRON STAR DETECTABILITY AND 
CONSTRAINTS 

To assess the detectability of the Galactic neutron star 
population, we use the methods outlined in Section |LT] to 
find the spin frequency v and gravitational wave strain 



4 Since we consider a continuous distribution in birth frequency, 
and a single star cannot be older than the Galaxy, neutron stars 
with high birth frequencies will not have existed long enough 
to have spun down past a certain frequency. Neutron stars will 
accumulate near this frequency causing a kink in the population's 
frequency distribution, as seen in [5]. 



amplitude h of every neutron star in our simulated pop- 
ulation. The population used in our analysis is described 
in detail in [5]. Although [5] presents three different ini- 
tial radial distributions, we choose to present our results 
using just the "gamma initial radial distribution" . We 
expect minimal deviation in our results if we were to 
consider the other two distributions presented in [Fj . The 
code that creates our simulated population provides the 
final position, which is easily turned into a distance from 
Earth r, and age t of each neutron star. Since a star's 
magnetic field strength B and ellipticity e dictate its fre- 
quency evolution, we must also choose what values of 
each to assign to the stars in our population. We fix e 
and B to a single value so that every neutron star in our 
population has the same e—B combination. While this 
approach will not mimic a realistic neutron star popula- 
tion, it is an important first step that provides valuable 
intuition for considering a more realistic population in 
the future. Lastly, we assign each star a birth frequency 
vq = 1/Poi where Pq is randomly drawn from the lognor- 
mal birth period distribution in [5]: 
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Here, Pq > 0.5 ms is the birth period, Pq = 5 ms is the 
mean, and a = 0.69 is the standard deviation. Given r, t, 
vq, e, and B, we can use the methods outlined in Section 
[Tljto find u(t, z^o | e, B) and h(r, v\e) for every neutron star 
in our simulated population. 

Once v and h are found for every neutron star, we 
can determine whether or not we expect a gravitational 
wave detector to detect our population. For simplicity, 
we only consider detection by a single Advanced Laser 
Interferometer Gravitational-wave Observatory (aLIGO) 
detector. We use the aLIGO noise curve for a single 
detector from [15] , which is the expected sensitivity of 
aLIGO as a function of gravitational wave frequency. To 
estimate the strain, we assume that we have a year of 
aLIGO data, and that the data is analyzed coherently in 
short 72 hr stretches, with the short stretches combined 
incoherently. This assumes the LIGO Scientific Collabo- 
ration will be doing similar searches to the ones currently 
done by Einstein@Home in the aLIGO era. An over- 
all trials factor of 100 is applied, which is considered a 
conservative estimate. We then compare each neutron 
star to the estimated noise curve to determine the num- 
ber of neutron stars in our population that aLIGO will 
be able to detect. We assume that a neutron star will be 
detected if its strain is above the aLIGO noise curve. To 
assess the detectability of the neutron star population, 
we construct the fraction 



(11) 



where N s i m is the number of stars in the simulated popu- 
lation, and ATjct is the number of stars aLIGO can detect 
from this population. To reduce statistical fluctuations, 
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FIG. 3: Contours of log 10 n (described in Section with respect 
to (log 10 B, log 10 e), which illustrate the expected detectability of 
the neutron star population with various e and B combinations. 
Our analysis was performed on populations with logarithmic spac- 
ings of 1/8 dex in e and B. The dashed lines are the analytic 
approximations for log 10 n described in Section |IV| We plot three 
results: log 10 n = —7,-8,-9, respectively from left to right. The 
dotted lines show the boundaries that separate where the analytic 
argument's assumptions are valid from where they are not (Section 
|IV| . They only hold for detectable neutron stars that are young 
(tmax 10 Myrs), found above the horizontal dotted line, and 
dominated by electromagnetic emission (r^ip < i~ gw ), found to the 
right of the positively sloped dotted line. 



we simulate many more neutron stars than are actually 
expected to be in our Galaxy. Multiplying this fraction n 
by the number of neutron stars in our Galaxy N sa i gives 
the number of detectable neutron stars in our Galaxy. If 
n-TVgai is greater than or equal to one, the population will 
likely be detectable; if it is less than one, the population 
will likely be undetectable. In Figure [3j we plot contours 
of log 10 n with respect to (log 10 B, log 10 e), illustrating 
the expected detectability of the neutron star population 
with various e and B combinations. 

We can further use our results (Figure [3]) to place 
bounds on the e-B parameter space of the Galactic neu- 
tron star population. In Figure [3j the contour corre- 
sponding to n ■ ./Vgai = 1 is the boundary above which 
the e-B parameter space is disallowed, assuming (pes- 
simistically) no aLIGO detection of continuous gravita- 
tional waves associated with a Galactic neutron star. In 
this way, Figure [3] sets lower bounds on B for fixed e val- 
ues (or upper bounds on e for fixed B values) if aLIGO 
does not make an isolated neutron star detection. For 
instance if N ga x ~ 10 9 [6], and we assume neutron stars 
have a typical ellipticity of e ~ 1CT 7 [6], Figure [3] shows 
that the minimum magnetic field strength of Galactic 
neutron stars is B > 10 11 G in the absence of an aLIGO 
detection. Conversely if N ga \ ~ 10 9 , and we assume 
neutron stars have a typical magnetic field strength of 
B ~ 10 11 G, Figure |3] shows that the population's maxi- 



FIG. 4: The top plot is a typical sensitivity curve of a ground-based 
gravitational wave detector. We assume that isolated neutron stars 
with /i(i'dct) > ^dct w iH be detected. Below are diagrams of the 
Milky Way disk. The black dot is Earth's location within the Milky 
Way disk. The two volumes V ma x outlined by dotted lines are the 
maximum volumes within which detectable neutron stars can be 
contained. Vmax will be the volume of a sphere if r max < 
(top diagram), and V ma x will be the volume of a spherical segment 
if »"max > Hmw (bottom diagram). See Section [TV] 



mum ellipticity is e < 10 -7 in the absence of an aLIGO 
detection. This argument also applies if aLIGO does 
make isolated neutron star detections. If 7V ga i ~ 10 9 , and 
we assume neutron stars have e ~ 10~ 7 , then the min- 
imum magnetic field strength of Galactic neutron stars 
is B > 10 10 G if aLIGO detects 10 neutron stars. Con- 
versely if iVg a i ~ 10 9 , and we assume neutron stars have 
B ~ 10 10 G, then the population's maximum ellipticity 
is e < 10~ 7 if aLIGO detects 10 neutron stars. 



IV. ANALYTIC RESULTS 

In all previous sections, we used numerical methods 
to assess the detectability of Galactic neutron stars and 
place constraints on the properties of the Galactic neu- 
tron star population. In this section, we present an ana- 
lytical approach to setting bounds on the e-B parameter 
space of the Galactic neutron star population. Bland- 
ford's analytic argument considers neutron stars that 
emit only gravitationally. Our analytic argument, while 
still simplistic, applies to stars dominated by electromag- 
netic emission. 

We first use the aLIGO sensitivity curve described 
in Section IHII to constrain the volume around Earth in 
which detectable neutron stars must be contained. If a 
neutron star is detected, it will tend to be at or near the 
detector's most sensitive frequency, which we call ^dct- 
For simplicity, we assume that a neutron star must have 
v ~ ^det to be detected. Therefore, a neutron star will be 
detected by a ground-based gravitational wave detector if 
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h(vdet) > ^dot, where hdet is the value of the strain ampli- 
tude for which the detector is most sensitive (see Figure 
Eh. This detectability condition, along with Equation |9j 
translates into a constraint on the distance from Earth 
of detectable neutron stars. The maximum distance r max 
at which a neutron star with v = u^et could be detected 
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The detectability condition also translates into a con- 
straint on the volume that encloses detectable neutron 
stars. First, we assume that neutron stars are born uni- 
formly throughout the Galactic stellar disk at a constant 
rate J\f, which is the number of births per unit time. The 
volume of the Milky Way, which we approximate to be a 
disk, is roughly 



Vmw = 7r^Mw(2-ffMw)> 



(13) 



where Rmw is the radius of the Galactic disk, and -Hmw 
is half its height (Figure |4| . The volume contained in 
fmax will be a sphere for r max < i?MW- However, for 



> Hi 



MW: 



the volume contained within 



will 



be a sphere with its top and bottom caps truncated by 
the top and bottom surfaces of the Milky Way disk, as 
illustrated in Figure [3] This shape is called a spherical 
segment. Therefore, the maximum volume K max in which 
neutron stars with v — v^e% could be detected is: 
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From the constraint on the volume that encloses de- 
tectable neutron stars of frequency v — Vdct-, we can find 
the minimum allowed magnetic field strength in the ab- 
sence of an isolated neutron star detection. Remember- 
ing our constant birth rate assumption, the average time 
i max between neutron star births into the volume V max is 
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assuming a uniform spatial distribution. In order to en- 
sure a neutron star detection, at least one star within the 
volume Vmax must have v > i>dct at all times. This will 
be the case if v{t max ) > Vdeti because a neutron star spin- 
ning down below Vdct will always be accompanied by a 
new star being born into the volume Vmax. Likewise, we 
also assume that, on average, when a detectable neutron 
star escapes Vmax due to its motion in the Galaxy, an- 
other detectable neutron star will enter Vmax- Assuming 
that a neutron star's frequency evolution is dominated 
by dipolar emission, we solve for the minimum magnetic 
field strength B min below which a neutron star detection 
is not guaranteed by substituting t ma x into Equation [4] 



and solving for B: 
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Neutron stars with B > S m in in V ma x will spin down to 
v < Vdct before another star is born into Vmax. There- 
fore, in the absence of an aLIGO detection, B = B m i n is 
the minimum possible magnetic field strength of Galactic 
neutron stars with fixed ellipticity values, since B < -B m in 
ensures a detection. 

Our argument is easily extended to the case of N s 
neutron star detections. To do this, we solve for when 

^(^s^max) > ^dct- The result adds a factor of N s 1 in 
front of Equation [T6| 

We have made several assumptions in setting up our 
analytical argument. It is important to emphasize two 
of our argument's most crucial assumptions in order to 
clearly outline the physical systems for which our argu- 
ment holds. The first crucial assumption is that the spa- 
tial distribution of neutron stars in the Milky Way is a 
uniform cylinder. Neutron stars will diffuse out of the 
Galactic disk due to Galactic acceleration and their kick 
velocities. The timescale for this process is found by di- 
viding the average kick velocity in [5] by the gravitational 
acceleration (found by dividing the gravitational poten- 
tial in [5] by the length scale). Therefore, our argument 
holds when t max < 10 Myrs. The second crucial assump- 
tion is that the frequency evolution of neutron stars is 
dominated by dipolar emission. Therefore, for Uq v, 
our argument holds when Tdi p i~ gw . 



I n Fi gure [3} we have plotted the relationship in Equa- 
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(with the factor of N s in front) as dashed lines 
on top of our numerical results. We use z/<jet ~ 100 Hz 
and hdet ~ 6.0 x 10~ 26 , which approximately corre- 
sponds to aLIGO 's most sensitivity strain and associ- 
ated frequency, and i?MW ~ 15 kpc, Hmw ~ 75 pc, and 
vq ~ 850 Hz, where -Rmw, Hmw, and Vq are estimated 
averages of the spatial and period distributions in [5] 
found by reducing the maximum values by a factor of e . 
We also use N ~ 0.02 years" 1 [14 . Our numerical re- 
sults should roughly follow these dashed lines, which cor- 
respond to n = 10~ 7 , 10~ 8 , 10~ 9 , respectively from left to 
right. The analytic results only hold for detectable neu- 
tron stars that are young (t maK < 10 Myrs), which cor- 
responds to the region above the horizontal dotted line, 
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FIG. 5: Extends our analytic argument to any gravitational wave 
detector. Equation |16| is plotted as solid contours of log 10 B with 
respect to (log 10 /idctjlogio e )i where B is in units of Gauss. The 
dashed and dotted lines show the boundaries that separate where 
our analytic argument's assumptions are valid from where they 
are not. Our argument does not hold in the gray, shaded re- 
gions; our argument does hold for neutron stars that are young 
(£max 10 Myrs), found below the dotted line (t max = 10 Myrs), 
and dominated by electromagnetic emission (r^p < T gw ), found 
below the dashed line (rdi p = Tgw)- 



and dominated by electromagnetic emission (T dip < r gw ), 
which corresponds to the region to the right of the pos- 
itively sloped dotted line. There is good agreement be- 
tween our numerical and analytic results, except in the 
transition region near the dotted boundaries where the 
analytic assumptions start to lose their validity. While 
the rough numerical values chosen for the parameters in 
our analytic argument can change the overall normaliza- 
tion of the analytic curves, the shape of the curves closely 
match the shape of the numerical contours. 

While we only consider detection by aLIGO in our nu- 
merical analysis, our analytical approach easily extends 
to any gravitational wave detector. Notice that in Equa- 
tion [16] £?min is a function of e, /idet, and Vdct- Therefore, 
we fix Vdet and plot contours of log 10 B nlin with respect 
to (log 10 e, log 10 ft-dct) in Figure [5] to illustrate how our 
argument extends to other detectors. 

It also seems natural to extend our argument to the 
gravitar case, in which the frequency evolution of neu- 
tron stars is dominated by gravitational emission, by 
solving Equation [2] for B m [ n under the assumption that 
ig W ^ Tdip. However, detectable gravitars can be older 
than 10 Myrs, thus violating our assumption that i max ^ 
10 Myrs. Therefore, these methods cannot be applied to 
the gravitar case. 



V. CONCLUSION 

We have used the methods described in Section [TT] to 
find the gravitational wave amplitude and spin frequency 
of every neutron star in the simulated population de- 
scribed in [5]. This involved allowing for both electro- 
magnetic and gravitational emission in a neutron star's 
frequency evolution (Equation [l]). We then solved for 
each neutron star's frequency using either Equation [2] 
[4j or [8] and each neutron star's gravitational wave strain 
amplitude using Equation [9j We used the simulated pop- 
ulation to assess the dctcctability of and set bounds on 
the e-B parameter space of the Galactic neutron star 
population. Our results are summarized in Figure [3] As- 
suming that the Galactic neutron star population con- 
sists of ./Vgai ~ 10 9 stars, and assuming aLIGO does not 
make a neutron star detection, the contour log 10 n = —9 
in Figure [3] separates the allowed e-B parameter space 
(below the contour) from the disallowed e-B parameter 
space (above the contour). In other words, assuming we 
know the magnetic field strength of the neutron star pop- 
ulation, we can place upper bounds on the population's 
ellipticity; or, assuming we know the ellipticity of the 
neutron star population, we can place lower bounds on 
the population's magnetic field strength. 

In this paper, we have only considered the simple (and 
unrealistic) case in which all neutron stars have the same 
magnetic field strength and ellipticity. However, we have 
demonstrated that both a gravitational wave detection or 
the lack of a gravitational wave detection can be used to 
constrain some of the properties of the Galactic neutron 
star population. To make confident quantitative state- 
ments regarding the properties of the Galactic neutron 
star population, we must construct a more realistic popu- 
lation. Moving forward, we plan to incorporate magnetic 
field and ellipticity distributions and evolutions into our 
analysis to more closely mimic the Galactic neutron star 
population [6][l5][l6]. 
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